Astronomy & Astrophysics manuscript no. He061 


February 5, 2008 


(DOI: will be inserted by hand later) 





The projection factor of 6 Cephei 

A calibration of the Baade-Wesselink method using the CHARA Array 

Antoine Merand 1 *, Pierre Kervella 1 , Vincent Coude du Foresto 1 , Stephen T. Ridgway 1 ' 2 ' 3 , Jason P. Aufdenberg 2 , 
Theo A. ten Brummelaar 3 , David H. Berger 3 , Judit Sturmann 3 , Laszlo Sturmann 3 , Nils H. Turner 3 , and Harold A. 

McAlister 3 

1 LESIA, UMR8109, Observatoire de Paris, 5, place Jules Janssen, 92195 Meudon (France) 

2 National Optical Astronomical Observatory 950 North Cherry Avenue, Tucson, AZ 85719 (USA) 

3 Center for High Angular Resolution Astronomy, Georgia State University, P.O. Box 3965, Atlanta, Georgia 30302-3965 
(USA) 

Received — / Accepted — 

Abstract. Cepheids play a key role in astronomy as standard candles for measuring intergalactic distances. Their distance is 
usually inferred from the Period-Luminosity relationship, calibrated using the semi-empirical Baade-Wesselink method. Using 
this method, the distance is known to a multiplicative factor, called the projection factor. Presently, this factor is computed using 
numerical models - it has hitherto never been measured directly. Based on our new interferometric measurements obtained with 
the CHARA Array and the already published parallax, we present a geometrical measurement of the projection factor of a 
Cepheid, S Cep. The value we determined, p = 1.27 ± 0.0 6, confirms the genera lly adopted value of p = 1.36 within 1.5 sigmas. 
Our value is in line with recent theoretical predictions of lNardetto et alj I2004h . 
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1. Introduction 

Cepheid stars are commonly used as cosmological distance 
indicators, thanks to their well-established Period-Luminosity 
law (P-L). This remarkable property has turned these super- 
giant stars into primary standard candles for extragalactic dis- 
tance estimations. With intrinsic brightnesses of up to 100,000 
times that of the Sun, Cepheids are easily distinguished in dis- 
tant galaxies (up to about 30 Mpc distant). As such, they are 
used to calibrate the secondary distance indicators (supernovae, 
etc..) that are used to estimate even larger cosmological dis- 
tances. For instance, the Hubble Key Pro ject to measure the 
Hubble constant Hq dFreedman et aljboOll) is based on the as- 
sumption of a distance to the LMC that was established primar- 
ily using Cepheids. Located at the very base of the cosmologi- 
cal distance ladder, a bias on the calibration of the Cepheid P-L 
relation would impact our whole perception of the scale of the 
Universe. 

1.1. Period-Luminosity calibration 

The P-L relation takes the form log L — a log P + /?, where L 
is the (absolute) luminosity, P the period, a the slope, and f3 
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the zero point. The determination of a is straightforward: one 
can consider a large number of Cepheids in the LMC, located 
at a common distance from us. Calibrating the zero-point /3 is 
a much more challenging task, as it requires an independent 
distance measurement to a number of Cepheids. Ideally, one 
should measure directly their geometrical parallaxes, in order 
to obtain their absolute luminosity. Knowing their variation pe- 
riod, j6 would then come out easily. However, Cepheids are rare 
stars: only a few of them are located in the solar neighborhood, 
and these nearby stars are generally too far away for precise 
parallax measurements, with the exception of 6 Cep. 

1.2. The Baade-Wesselink method 

The most commonly used alternative to measure the distance 
to a pulsating star is the Baade-Wesselink ( BW) method . 
Developed in th e first part of the 20th century teaaddll92d 
IWesseHnkll^ . it utilizes the pulsational velocity V pu ] s of the 
surface of the star and its angular size. Integrating the pulsa- 
tional velocity curve provides an estimation of the linear radius 
variation over the pulsation. Comparing the linear and angular 
amplitudes of the Cepheid pulsati on gives directly its d istance. 
The most recent implementation jKervella et alj|2004l) of the 
BW method makes use of long -baseline interferometry to mea- 
sure directly the angular size of the star. 
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Unfortunately, spectroscopy measures the apparent radial 
velocity V ia ±, i.e. the Doppler shift of absorption lines in the 
stellar atmosphere, projected along the line of sight and inte- 
grated over the stellar disk. This is where p, a projection factor, 
has to be introduced, which is defined as p — Vpuis./Kad.- The 
general BW method can be summarized in the relation: 

XT 
(V md .(t)-Vy)dt (1) 

where d is the distance, p the projection factor, 6 the angu- 
lar diameter and V y the systematic radial velocity. There are in 
fact many contributors to the p-factor. The main ones are the 
sphericity of the star (purely geometrical) and its limb darken- 
ing (due to the stellar atmosphere structure). A careful theo- 
retical calculation of p requires modeling dynamically the for- 
mation of th e absorption line in the pulsating atmosphere of 
the C epheid dParsonsll 1 9721 Is abbey et alJll995UNardetto et alJ 
2004). 

Until now, distance measurements to Cepheids used a p- 
factor value estimated from numerical models. Looking closely 
at Eq. H it is clear that any uncertainty on the value of p will 
create the same relative uncertainty on the distance estimation, 
and subsequently to the P-L relation calibration. In other words, 
the Cepheid distance scale relies implicitly on numerical mod- 
els of these stars. But how good are the models? To answer this 
question, one should confront their predictions to measurable 
quantities. Until now, this comparison was impossible due to 
the difficulty to constrain the two variables 9{T) and d from 
observations, i.e. the angular diameter and the distance. 

Among classical Cepheids, 5 Cep (HR 8571, HD 213306) 
is remarkable: it is not only the prototype of its kind, but also 
the Cepheid with the most precise trigonometric parallax cur- 
rently available, obtained recently using the FGS instru ment 
aboard the Hubble Space Telescope Benedict et alJ2002l) . This 
direct measurement of the distance opens the way to the direct 
measurement (with the smallest sensitivity to stellar models) 
of the p factor of 6 Cep, provided that high-precision angular 
diameters can be measured by interferometry. 

2. Application of the BW method to 5 Cep. 

To achieve this goal, interferometri c observations were 
under t aken at the CHARA Arr ay Jten Brummelaar et alJ 
120031 Iten Brummelaar et al.l 120051) . in the infrared K' band 
( 1 .95 pm < A < 2.3 um) with the Fiber Linked U nit for Optical 
Recombination dCoude du Foresto et al.l Eo03l) (FLUOR) us- 
ing two East-West baselines of the CHARA Array: El-Wl 
and E2-W1, with baselines of 313 and 251 m respectively. 
Observations took place during summer 2004 for E2-W1 
(seven nights between JD 2 453 216 and JD 2 453 233) and Fall 
2004 for El-Wl (six consecutive nights, from JD 2 453 280 
to JD 2 453 285). The pulsation phase wa s computed using 
the fo llowing period and reference epoch JMoffett & Barnes! 
Il985l) : P = 5.3663 16 d, r = 2 453 674.144 (Julian date), the 
0-phase being defined at maximum light in the V band. The 
resulting phase coverage is very good for the longest baseline 
(El-Wl), while data lack at minimum diameter for the smaller 
one(E2-Wl) 





S. Type 


UD diam. (mas) 


Baseline 


HD 2952 


K0III 


0.938 ±0.013 


Wl-El 


HD 138852 


K0III-IV 


0.952 ±0.012 


Wl-El 


HD 139778 


Kim: 


1.072 ±0.014 


W1-E2 


HD 186815 


K2III 


0.713 ±0.009 


W1-E2 


HD 206349 


K1II-III 


0.869 ±0.011 


Wl-El, W1-E2 


HD 206842 


Kim 


1.214 ±0.016 


W1-E2 


HD 214995 


K0III: 


0.947 ±0.013 


Wl-El 


HD 216646 


K0III 


1.051 ±0.015 


Wl-El, W1-E2 


HD 217673 


K1.5II 


1.411 ±0.020 


W1-E2 



Table 1. Calibrators with spectral type, uniform disk an- 
gular diameter in K band (in milliarcsecond) and baseline 
jMerand et all2005l) . 



The FLUOR D ata reduction software (DRS) 

dCoude du Foresto et alJ Il997l) . was used to extract the 
squared modulus of the coherence factor between the two 
independent apertures. All calibrator stars wer e chosen in a 
catalo gue computed for this specific purpose JMerand et alJ 
2005) (see Table 0. Calibrators chosen for this work are all 
K giants, whereas 6 Cep is a GO supergiant. The spectral type 
difference is properly taken into account in the reduction, even 
though it has no significant influence on the final result. The in- 
terferometric transfer function of the instrument was estimated 
by observing calibrators before and after each 6 Cep data point. 
The efficiency of CHARA/FLUOR was consistent between 
all calibrators and stable over the night around 85%. Data 
that share a calibrator are affected by a common systematic 
error due to the uncertainty of the a priori angular diameter 
of this calibrator. In order to interpret our data properly, we 
used a specific formalism JPerrinll2003l) tailored to propagate 
these correlations into the model fitting process. Diameters are 
derived from the visibility data points using a full model of the 
FLUOR instrument including the spectral bandwidth effects 
jKervella et alJE()03l) . The stellar center-to-limb darkening is 
correct ed using a mo del intensity profile taken from tabulated 
values aretll2000D with parameters corresponding to <5 Cep 
(T e ff = 6000/ir, logg = 2.0 and solar metallicity). The limb 
darkened (LD) angular diameter comes out 3% larger than its 
uniform disk (UD) counterpart. 

The theoretical correction for LD has only a weak influ- 
ence on the p-factor determination, since that determination is 
related to a diameter variation. For example, based on our data 
set, a general bias of 5% in the diameters (due to a wrongly 
estimated limb darkening) leads to a bias smaller than 1% in 
terms of the p-factor. Differential variations of the LD correc- 
tion during the pulsation may also influence the projection fac- 
tor: co mparison between hyd rodynamic and hydrostatic simu- 
lations jMarengo et all2003l) showed negligible variations. An 
accuracy of 0.2% on the angular diameters for a given baseline 
is required to be sensitive to dynamical LD effects. This is close 
to, but still beyond, the best accuracy that we obtained on the 
angular diameter with a single visibility measurement: 0.35% 
(median 0.45%). 

Among the various sets of measurements of the radial ve- 
locity V,-ad.(0 available for 6 Cep, we chose measurements from 
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Fig. 1. Radial Velocity smoothed using splines. A. Radial ve- 
locity data points, as a function of pulsation phase (0-phase de- 
fined as the maximum of li ght). This set was ex tracted using a 
cross-correlation technique dBersier et all 19941) . The solid line 
is a 4-knot periodic cubic spline fit. B. Residuals of the fit. 



iBersier et all i ll 994 and lBarnes et all d2005l) . These works of- 
fer the best phase coverage, especially near the extrema, in or- 
der to accurately estimate the associated photospheric ampli- 
tude. In order not to introduce any bias due to a possible mis- 
match in the radial velocity zero-point between the two data 
sets, we decided to reduce them separately and then combine 
the resulting p-factor. An integration over time is required to 
obtain the photospheric displacement (see Eq[Q. This process 
is noisy for unequally spaced data points: the radial velocity 
profile was smoothly interpolated using a periodic cubic spline 
function. 

Fitting the inferred photospheric displacement and ob- 
served angular diameter variations, we adjust three parame- 
ters: the mean angular diameter 9, a free phase shift 0o and 
the projection factor p (see Fig. 1). The mean angular diam- 
eter is found to be 1.475 ± 0.004 mas (milliarcsecond) for 
both ra dial velocity data set s. Assuming a distance of 274 ± 
11 pc jBenedict et aljEoO^ . this leads to a linear radius of 
43.3 ±1.7 solar radii. The fitted phase shift is very small in 
both cases (of the order of .01). We used the same parame- 
ters I Moffett & Barnesll 19851) to compute the phase from both 
observation sets and considering that they were obtained more 
than ten years apart, this phase shift corresponds to an uncer- 
tainty in the period of approximately five seconds. We thus con- 
sider the phase shift to be reasonably the result of uncertainty 
in the ephemeris. 

The two different radial velocity data sets lead to a consoli- 
dated value of p — 1 .27 ± 0.06, once again assuming a distance 
of 274 ± 1 1 pc. The final reduced^ 2 is 1.5. The error bars ac- 
count for three independent contributions: uncertainties in the 
radial velocities, the angular diameters and the distance. The 
first was estimated using a bootstrap approach, while the oth- 
ers were estimated analytically (taking into account calibration 
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Phase 

Fig. 2. p-factor determination. A. Our angular diameter mea- 
surements (points). Crosses correspond to the medium base- 
line (E2-W1), while circles correspond to the largest base- 
line (El-Wl). The continuous line is the integration of the 4- 
knots periodic cubic spline fitted to the radial velocities (FigE) 
. Integration parameters: 9 = 1.475 mas, p = 1.269 and 
d = 21 A pc. B. Residuals of the fit 



p 


± C"Vrad. 


1.269 
1.280 


± 0.008 
± 0.012 


ref. (1) 
ref. (2) 


p 


± 0"Vrad. 


1.273 


± 0.007 


consolidated 




0"interf. 




± 0.020 






C dist. 




± 0.050 






p 


1.27 


± 0.06 





Table 2. Best fit results for p, with the two different radial 
velocity sets. The third line is a weighted average of the two 
individual measurements. Fourth and fith lines are the detailed 
quadratic contribution to the final error bar. Last line gives the 
fin al adopted value wit h the o verall error bar. Refe rences are: 
d) lBersieret"al](ll994 and (2) lBarnes etai1d2005h 



correlation for interferometric errors): for p, the detailed error 
is p = 1.273 ± 0.007 vrad. ± 0.020 interf . ± 0.050 di st.. The error is 
dominated by the distance contribution (see Tabled. 

3. Discussion 

Until now, the p-fa ctor has been dete rmined using models: 
hydrostatic models dBurki et alJ Il982l) produced the gener- 
ally adopted value, p - 1.36. First attempts were made by 
ISabbev et al.ldl995l) to take into account dynamical effects due 
to the pulsation. They concluded that the average value of p 
should be 5% larger than in previous works (1.43 instead of 
1.36) and that p is not constant during the pulsation. Because 
they increased p by 5%, they claimed that distances and diam- 
eters have to be large r in the same proportion. More recently 
iNardetto et alJ d2004 computed p specifically for S Cep using 
dynamical models. Different values of p were found, whether 
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one measures diameters in the continuum or in the layer where 
the specific line is formed. In our case, broad band stellar inter- 
ferometry (angular diameters are measured in the continuum) 
these authors suggest p = 1.27 ± 0.01. Concerning the vari- 
ation of p during the pulsation, they estimate that the error in 
terms of distance is of the order of 0.2%, smaller than what 
we would have been able to measure with our interferomet- 
ric data set. While our estimate, p = 1.27 ± 0.06, is statis- 
tically compatible with this recent work, marginally with the 
widely used p = 1.36, and not consistent with the former value 
p = 1.43 at a 2 cr level. We note that Gieren et al. (2005) have 
recently derived an expression of the p-factor as a function of 
the period that predicts a value of 1.47 + 0.06 for 5 Cep. While 
this va lue is in agreement with the modeling by ISabbev et alJ 
l l 19951) . is is slightly larger than the present measurement (by 
2.4 cr). As a remark, Gieren et al. obtain a distance of 280±4pc 
for 5 Cep, that is slightly larger than Benedict et al.'s (2002) 
value 274 + 1 1 pc assumed in the present work. Assuming this 
new distance estimation with our data would result in a ^-factor 
of 1.30 + 0.06, bringing the agreement to 2 cr only. 

Our geometrical determination of the p-factor, p — 1 .27 ± 
0.06, using the IB W method is currently limited by the error 
bar on the parallax jBenedict et aljE002i) . Conversely, assum- 
ing a perfectly known p-factor, the uncertainty of the stellar 
distance determined using the same method would have been 
only 1.5%, two-times better than the best geometrical parallax 
currently available. The value we determined for p is statisti- 
cally compatible with the value generally adopted to calibrate 
the Cepheid P-L relation in most recent works. It is expected 
that the distance to approximatively 30 Cepheids will be deter- 
mined interferometrically in the near future usin g particularly 
the C HARA Array and the VLT Interferometer ( iGlindemannl 
120051) . In order not to limit the final accuracy on the derived 
distances, theoretical p-factor studies using realistic hydrody- 
namical codes is necessary. With a better understanding of the 
detailed dynamics of the Cepheid atmospheres, we will be in a 
position to exclude a /^-factor bias on the calibration of the P-L 
relation, at a few percent level. 
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Table 3: Individal measurements, colums are (1) date of observation, 
JDo=2 453 200.5 (2) phase (3,4) u-v coordinate in meter (5) squared vis- 
ibility and error (6) corresponding limb darkened disk diameter in mas 
(7,10) HD number of calibrators, prior and after the given data point re- 
spectivaly, means that there was no calibrator (8 ,9,11,12) quantities 
for computing the correlation matrix llPerrinll2003h : cr v i are errors on 
the estimated visibility of the calibrators. NOTE : this table will be 
archived electronically 



